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This paper presents cosmological results from the final data release of the WiggleZ Dark Energy 
Survey. We perform full analyses of different cosmological models using the WiggleZ power spectra 
measured at z=0.22, 0.41, 0.60, and 0.78, combined with other cosmological datasets. The limiting 
factor in this analysis is the theoretical modelling of the galaxy power spectrum, including non- 
linearities, galaxy bias, and redshift-space distortions. In this paper we assess several different 
methods for modelling the theoretical power spectrum, testing them against the Gigaparsec WiggleZ 
simulations (GiggleZ). We fit for a base set of 6 cosmological parameters, {Q,hh? , JIcdm/i^, ^^o, t, 
Ag, rig}, and 5 supplementary parameters {nmn, r, w, flk, "^m^}. In combination with the Cosmic 
Microwave Background (CMB), our results are consistent with the ACDM concordance cosmology, 
with a measurement of the matter density of fim = 0.29 ± 0.016 and amplitude of fiuctuations as ~ 
0.825 ± 0.017. Using WiggleZ data with CMB and other distance and matter power spectra data, 
we find no evidence for any of the extension parameters being inconsistent with their ACDM model 
values. The power spectra data and theoretical modelling tools are available for use as a module 
for CosmoMC, which we here make publicly available at http://smp.uq.edu.au/wigglez-datal 
We also release the data and random catalogues used to construct the baryon acoustic oscillation 
correlation function. 



PACS numbers: 98.80.-k 

I. INTRODUCTION 

The aim of the WiggleZ Dark Energy Survey has been 
to measure the large scale structure of the Universe and 
use that to learn about the properties of dark energy. 
Here we present our final data release and cosmology re- 
sults, using the full large scale structure power spectrum. 
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The Universe is filled with structure, forming under 
gravitational collapse. We observe this structure through 
light, emitted when the gas is compressed and heated. 
Cosmology is the study of this light, from stars in dis- 
tant galaxies, to cosmic microwaves emitted by the Big 
Bang - and inferences are made as to the structure that 
it traces. We have built up a complex picture of struc- 
ture formation, where it is seeded as tiny density fluc- 
tuations and grows under gravitational collapse to form 
the 'cosmic web' of galaxies and clusters that we observe 
in the late-time Universe. The cosmological model of 
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acausal-seeded, hierarchical structure formation in a ge- 
ometrically flat, late-time accelerating universe has been 
well tested using precision cosmological datasets such 
as measurements of the Cosmic Microwave Background 
(CMB) by the COsmic Background Explorer (COBE) 
and Wilkinson Microwave Anisotropy Probe (WMAP), 
the expansion history of the Universe using standard can- 
dles and rulers such as Type-la Supernovae (SNIa) and 
Baryon Acoustic Oscillations (BAO), and the late Uni- 
verse matter power spectra (P(fc)) and growth of struc- 
ture. The name of this 'concordance' cosmological model 
is ACDM. 

The definition of a cosmological model is one that de- 
scribes the whole Universe, and possesses a number of 
tuneable parameters. In the past, cosmological models 
were incommensurable, in that it was almost impossible 
to judge one in the light of another (such as the Big Bang 
model compared to the Steady State model). The ad- 
vent of precision cosmological data has allowed ACDM 
to emerge as the predominant model, but it must al- 
ways be checked to be consistent with new data as it be- 
comes available. New alternatives/extensions must also 
be tested, to see if they make better predictions than the 
standard model, which would indicate that it should be 
updated. 

The concordance cosmological model is described by 
six free parameters: the Hubble parameter today, Hq 
(where the dimensionless h = i?o/100kms~^Mpc~^), the 
physical baryon density, flhh'^, the physical Cold Dark 
Matter density, ilcDM^i^, the optical depth to reionisa- 
tion, r, the amplitude of the primordial scalar density 
perturbations on scales oi k — 0.05/i"^Mpc, Ag, and 
the spectral index of the primordial power spectrum of 
density perturbations, Ug. Extensions to this standard 
model involve adding one or more extra parameters such 
as curvature fi^, the equation of state of the dark en- 
ergy, w, the mass of the neutrinos, X]"^''' running 
of the spectral index of the primordial power spectrum 
of density perturbations, riiun (= dlnn^/dlnfe), and the 
amplitude of the primordial tensor perturbations (gener- 
ated by primordial gravitational waves from the end of 
inflation), r. All of these represent some extra physical 
effects that have been suggested theoretically and would 
contribute to the expansion and/or structure formation 
of the Universe, but may not be detectable using current 
cosmological data. 

The advantage of a complete model of structure forma- 
tion is that it allows us to combine information from the 
early and late Universe to test cosmological models. This 
has been done by the CfA Redshift Survey [1 , the APM 
Galaxy Survey 0, VIMOS-VLT Deep Survey (VVDS) 
[3], the Two-Degree Field Galaxy Redshift Survey (2dF- 
GRS) 15 and the Sloan Digital Sky Survey (SDSS) |5]. 

In this paper we present analyses of the cosmological 
models using data from the WiggleZ Dark Energy Survey 
[6]. WiggleZ is a large-scale galaxy redshift survey con- 
ducted with the AAOmega multi-object spectrograph on 
the Anglo- Australian Telescope at Siding Springs Obser- 



vatory. WiggleZ was designed to study the effect of dark 
energy on the expansion history of the Universe and on 
the growth of cosmological structures across an unprece- 
dented period of cosmic history. We have already pre- 
sented our measurements of several important features 
such as the baryon acoustic oscillation scale [71 [8j, the 
growth of structure 9], the Alcock-Paczynski effect, [lOj . 
and the mass of the neutrino llj. 

An important step in these types of analysis is generat- 
ing accurate predictions from theory for what the large- 
scale structure in the Universe should look like. In prac- 
tice that usually means predicting the shape of the galaxy 
power spectrum or correlation function. Using only linear 
theory, we can easily make predictions using perturbation 
theory. However, the relationship between the observed 
galaxy power spectrum and the linear perturbation the- 
ory prediction for the matter power spectrum is com- 
plicated by three physical effects: non-linear structure 
formation, galaxy bias, and redshift space-distortions. 

Non-linear structure formation: When density pertur- 
bations become large enough they are no longer well de- 
scribed by flrst order perturbation theory. For testing 
cosmological models using the CMB power spectrum, the 
linear theory prediction is accurate enough when com- 
pared to the data (at least at the moment), as the am- 
plitude of the density fluctuations at z ~ 1100 are very 
small. For the matter power spectrum at low-redshift 
this is not the case, as some non-linear evolution of the 
density field has taken place. The amplitude of this non- 
linear evolution can be estimated by numerical simula- 
tions [I2HI4] . and then applied as a correction to the 
linear prediction (e.g., Halofit [T5]l. 

Galaxy bias: This is not the end of the story, since 
we observe the distribution of galaxies, but our the- 
ory predicts the distribution of matter, which includes 
non-luminous, non-baryonic dark matter. The galax- 
ies are slightly decoupled from the matter, through the 
complex baryonic physics of star and galaxy formation. 
The simplest coupling to be assumed is the idea of the 
linear 'bias' [16] . an overall, shape- independent ampli- 
tude scaling from the matter power spectrum to the 
galaxy power spectrum. This bias parameter is differ- 
ent for different populations of galaxies, and relates di- 
rectly to the galaxy-formation history of the population 
being sampled. Because of this, we also expect the bias 
to evolve with redshift and environment (e.g., density 
where the galaxies are found), which would introduce 
scale-dependence to the bias. Using numerical simula- 
tions of galaxies, it should be possible to predict the bias 
for a given population of galaxies, and the technology has 
advanced to such a state that this is becoming possible. 

Redshift space distortions: Another complication arises 
because all our observed distances are inferred from mea- 
surements of redshifts (and angular separations). How- 
ever, galaxy redshifts do not solely come from the ex- 
pansion of the Universe, but also local velocities induced 
by gravitational over-densities. These peculiar velocities 
introduce redshift-space distortions into the clustering 
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power spectrum (as measured using the WiggleZ data in 
[H [To] ) , and will have a measurable effect on scales small 
enough for peculiar velocities to be significant. The effect 
of redshift space distortions can be detected by observing 
differences in the clustering pattern measured perpendic- 
ular to and parallel to the line-of-sight. 

A key part of the analysis in this paper is a detailed 
exploration of the different approaches we have used to 
successfully model the power spectrum in terms of the 
cosmological parameters we are constraining. We have 
tested seven different approaches for our analysis. The 
models are fairly similar at low values of the comoving 
wavenumber k (large scales), where the large-scale clus- 
tering can be treated as linear. There the theory is quite 
robust, and there is little difference between the predic- 
tions. However the difference between the models starts 
to increase for k > 0.2 /i Mpc~^, which significantly af- 
fects the outcome of the fitting. As we discuss in section 



IV D| we are fitting up to fcmax — 0.3 hMpc~ , as we 
found this to be the optimal compromise between robust 
modelling and data inclusion. 

With this paper, we also make public our data, as well 
as a code to implement the preferred power spectrum 
modelling method. The power spectrum data (measure- 
ments, covariance matrices and window functions) we 
make available along with the modelling code as a module 
for the CosmoMC cosmological analysis package [17] . We 
also make available the data that was used to construct 
the baryon acoustic peak correlation function measure- 
ments reported in Blake et al. [S]. 

The structure of the paper is as follows. In Section 
jn] we describe the WiggleZ Dark Energy Survey, the 
data that is available, the selection functions, the A^-body 
simulations (GiggleZ), and the random catalogues — in 
other words, all the raw material you need in order to 
perform the analysis in this and our other papers. We 
then describe in Section IIIII how we extracted the ob- 
served power spectrum from this data. In Section |IV| 
we describe how we model the galaxy power spectrum 
theoretically, starting with the linear matter power spec- 
trum and making additions to account for extra physical 
effects. In particular in this section we discuss the dif- 
ferent formalisms for dealing with non-linear structure 
formation, galaxy bias, and redshift space distortions, 
and select a particular model to perform our analyses. 
In Section |V] we conduct statistical analyses of a number 
of different cosmological models, combining the WiggleZ 
data with independent observations of the CMB, BAO, 
LSS, and SNIa. Finally, we summarise our results in 
Section 



II. THE WIGGLEZ DARK ENERGY SURVEY 
A. Data 

The WiggleZ Dark Energy Survey was designed to de- 
tect the BAO scale at higher redshifts than was possible 
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FIG. 1: (color online). The WiggleZ power spectrum in the 
four redshift bins, weighted averaged over the seven regions. 
We artificially offset the redshift bins in the vertical direc- 
tion to distinguish them (with the lowest redshift bin at the 
bottom moving up to highest at the top). At each redshift 
the line gives the best fit flat-ACDM cosmology prediction, 
convolved with the window function for each region and then 
averaged over the regions. The shaded region gives the range 
in wavenumber of data we used for the analysis. 



with previous datasets. Galaxies were initially selected 
by a joint selection from optical galaxy surveys (from 
the Sloan Digital Sky Survey [5] in the Northern Galac- 
tic Cap, and from the Red Sequence Cluster Survey 2 
(RCS2) [IH] in the Southern Galactic Cap) and from ul- 
traviolet imaging by the Galaxy Evolution Explorer satel- 
lite (GALEX). A number of magnitude and colour cuts 
were made to select high-redshift star-forming galaxies 
with prominent emission lines, and then spectra of these 
galaxies were taken by the AAOmega spectrograph [12] 
in 1-hour exposures at the Anglo- Australian Telescope. 
In total 238, 000 galaxy redshifts were measured by the 
survey. 

These galaxies were used to map seven regions of the 
sky with a total volume of 1 Gpc'^ in the redshift range 
z < 1 [6]. We split the data into four redshift bins with 
ranges 0.1 < z < 0.3, 0.3 < z < 0.5, 0.5 < z < 0.7 
and 0.7 < z < 0.9. The number of galaxies included in 
each redshift slice is listed in Table |lj We determined the 
effective redshift of the power spectrum estimate for each 
slice by weighting each pixel in the selection function by 
its contribution to the power spectrum error, as given by 
Eq. (13) in [7 . The effective redshifts of these bins is 
evaluated to be z^s = [0.22,0.41,0.6,0.78]. The power 
spectra, Pobsj and covariance matrices, C, are measured 
in Ak = 0.01 /iMpc"^ bins using the optimal-weighting 
scheme proposed by Feldman et al. j20] for a fiducial 
cosmological model with matter density = 0.27 |21| . 
These power spectra and covariance matrices, along with 
the window functions for the power spectrum, we make 
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TABLE I: Number of galaxies used in the power spectrum 
estimation, divided by survey region and redshift slice. 



Region 0.1 < 2 < 0.3 0.3 < 2 < 0.5 0.5 < z < 0.7 0.7 < z < 0.9 



9-lir 


2671 


7438 


11294 


5189 


11-hr 


3446 


8118 


11766 


5910 


15- hr 


3578 


10942 


17928 


8943 


22-hr 


5782 


7728 


8924 


5337 


0-hr 


2814 


2986 


5667 


3689 


1-hr 


2654 


2596 


5016 


3144 


3-hr 


2973 


3522 


5982 


4315 



available as part of our cosmology analysis package.^ 

The region-averaged power spectrum for each of the 
redshift bins, and the best fit fiat ACDM cosmology curve 
with these window function effects taken into account, 
are shown in Fig. [T] and the covariance matrices in Fig. 

m 

The methodology used to construct the WiggleZ sur- 
vey selection function and random catalogues is described 
in Blake et al. [5T]. Monte Carlo realizations of cata- 
logues are produced following the angular and redshift 
dependences of the selection function for each region. 
This process models several effects including the vari- 
ation of the GALEX target density with dust and ex- 
posure time, the incompleteness in the current redshift 
catalogue, the variation of that incompleteness imprinted 
across each 2-degree field by constraints on the position- 
ing of fibres and throughput variations with fibre posi- 
tion, and the dependence of the galaxy redshift distribu- 
tion on optical magnitude. 



1. BAO Data 

In conjunction with this paper, we make available 
the data that we used to construct the baryon acous- 
tic peak correlation function measurements reported in 
Blake et al. [5]. This dataset consists of measurements 
in three overlapping redshift slices 0.2 < z < 0.6, 
0.4 < z < 0.8 and 0.6 < z < 1.0 for the 9-hr, 11-hr, 15- 
hr, 22-hr, 1-hr and 3-hr survey regions. For each region 
and redshift slice we make available the data catalogue 
and 10 random catalogues. We also provide the result- 
ing correlation function and covariance matrix for each 
measurement (derived from lognormal realizations) . The 
number of galaxies included in each data subset is listed 
in Table [nj (Note that in the BAO measurement we do 
not use the 0-hr region, due to incompleteness concerns). 



The package contains the data combined with 
a module for CosmoMC, and is available at 
http: // smp. uq. edu. au/wigglez-data. 



TABLE II: Number of galaxies divided by survey region and 
redshift slice used in the baryon acoustic peak correlation 
function analysis reported by Blake et al. [8|. 



Region 0.2 < z < 0.6 0.4 < z < 0.8 0.6 < z < 1.0 



9-hr 


15128 


18978 


11424 


11-hr 


19202 


23940 


15064 


15-hr 


22309 


30015 


19471 


22-hr 


15884 


16146 


11024 


1-hr 


6927 


9437 


7880 


3-hr 


8000 


10241 


8756 



This sample is constructed from the full WiggleZ red- 
shift catalogue in each of these regions by applying cuts 
for redshift and contiguity. 

Note that this 'Large-Scale Structure catalogue' is a 
sub-set of the full redshift catalogue, but no matching 
random catalogues exist for the full redshift catalogue, 
only for this LSS sample. This only covers the sample 
used for the BAO measurement, and cannot be used for 
computing the full galaxy power spectrum and its errors. 



B. The GiggleZ simulation 

The Gigaparsec WiggleZ Survey (GiggleZ) simulations 
[22] were designed to probe the low-mass haloes traced 
by WiggleZ galaxies, whilst providing an equivalent sur- 
vey volume allowing the measurement of power spectrum 
modes with k = 0.01- 0.5 /iMpc~^. They provide a 
powerful means for testing and calibrating our model- 
ing algorithms. The main simulation we use is a 2160^ 
particle dark matter simulation run in a periodic box 
of side Ih^^ Gpc. The resulting particle mass of this 
simulation is 7.5 x 10^h~^ Mq which permits us to re- 
solve bound systems with masses ~ 1.5 x lO^^h^^ Mq, 
facilitating studies of halos with clustering biases near 
unity such as WiggleZ galaxies. The initial conditions of 
the simulation were constructed using a CAMB^ power 
spectrum based on a flat ACDM cosmology using the 
best fit WMAP 5-year values with {$7^, ^^b, /i, erg, nj} = 
{0.273, 0.0456, 0.705, 0.812, 0.960}. 

Bound structures in the dark matter catalogue were 
identified using Subfind [33], which uses a friends-of- 
friends (FoF) scheme followed by a substructure analy- 
sis to identify bound overdensities within each FoF halo. 
We use the Subfind substructures for our analysis, tak- 
ing the value of each halo's maximum circular velocity 
as a proxy for mass. In order to generate mock WiggleZ 
catalogues at a particular redshift we selected the range 
of halo groupings with large-scale clustering bias closest 
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FIG. 2: (color online). The WiggleZ covariance matrix in the 

to the WiggleZ sample under analysis, and applied the 
WiggleZ selection function to the mock dataset. 

By comparing simulated and observed power spec- 
tra, we found that over the range of scales and halo 
masses relevant for this analysis, the galaxy bias is scale- 
independent to within 1% [22] . 



III. TAKING FIDUCIAL MODEL AND 
WINDOW FUNCTION INTO ACCOUNT 

The method we use for generating a power spectrum 
from our measured positions of galaxies has been detailed 
in Blake et al. ^21,, and here we expand it to the entire 
WiggleZ sample. Crucially, the observed power spectrum 
is naturally convolved with the survey window function. 
(The window function quantifies the holes and incom- 
pleteness in the data that are due to the observing strat- 
egy rather than any real large scale structure.) In ad- 
dition, when generating this observed power spectrum 
we have used a fiducial cosmological model to convert 
observed angles and redshifts into distances. Before the 
observed galaxy power spectrum can be compared to the- 
oretical models, we must take those factors into account. 

The fiducial model and window function are accounted 




0.1 0.2 0.3 0.4 



^2 

IT redshift bins, as a weighted average over the seven regions. 

for by convolving the theoretical power spectrum with 
the survey window function, and scaling all distances 
with respect to the fiducial model 

Details of the window function, Wij{k), can be found in 
section [TT| and Blake et al. [2Tj. The scaling, Ogci, takes 
into account that the observed galaxy redshift-space posi- 
tions are converted to real space positions using a fiducial 
cosmology. It is calculated as [24l - f26] 

[DaMzosWH{z,s)' ^ > 

where Da is the angular diameter distance and H is 
the Hubble parameter evaluated for the effective redshift, 
Zoff, of the galaxy sample. A subscript 'fid' means this is 
the value in the fiducial model, while no subscript means 
this is the value in the model being tested. For a ACDM 
Universe Da{z) is given by 

(3) 
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where S is sin(a;), x or sinh(a;) depending on whether the 
Universe is closed, flat or open respectively, and H{z) 
given by the Friedman equation, 

^^-Jfi = a-,o(l + 2)^ + a„,o(l + z)' + f^fc,o(l + 2)2 + r!A,o , 
-"o 

(4) 

which can be integrated numerically. 



IV. MODELLING THE GALAXY POWER 
SPECTRUM 

As mentioned in the introduction, the limiting fac- 
tor in extracting cosmological constraints from a galaxy 
power spectrum lies in accurately modelling the theo- 
retical power spectrum in the cosmological model being 
tested. A major part of the analysis of this paper has 
been testing several different approaches to modelling the 
power spectrum. In testing this modelling we relied heav- 
ily on comparison to iV-body simulations (GiggleZ) to 
validate our treatment of non-linear structure formation, 
galaxy bias, and redshift space distortions. 

In terms of these modelling challenges the WiggleZ 
Dark Energy Survey has several potential advantages 
over previous surveys: 

Higher redshift: The effect of non-linear structure 
formation increases with time as density perturbations 
grow. For the distant galaxies probed by WiggleZ, the 
contamination from non-linearities is smaller than for 
previous surveys. 

Lower scale-dependence in bias: The relationship 
(bias) between the observed galaxy distribution and the 
dark matter distribution depends on the observed galaxy 
type. Previous studies [e.g. [171 [2S] measured red galaxies, 
which tend to cluster in the center of dark matter halos, 
whereas the star-forming blue WiggleZ galaxies avoid the 
densest regions. This leads to a lower overall bias, which 
makes WiggleZ less susceptible to any possible systemat- 
ics that could arise from a scale-dependence of the bias. 

Low non-linear pairwise velocities: The WiggleZ galax- 
ies have less significant fingers-of-god effects (due to lower 
non-linear pairwise velocities) because these low-mass 
tracers avoid massive clusters. 



A. Seven approaches to modelling P{k) 

Here we describe seven different approaches to mod- 
elling the power spectrum. They differ primarily in the 
way they treat non-linearities and bias. In the next sec- 
tion we will test these approaches by how well they re- 
produce the input cosmology of the GiggleZ simulation. 



1. Model A: Linear 

We use CAMB to calculate the linear matter power 
spectrum Pm{k) at redshift z for each set of cosmological 
parameters. We assume a linear bias b given by: 

Pg.lik)=b^Prn{k). (5) 

Since the linear bias s constant with respect to k, we 
analytically marginalise over b [17], 

Epj /~<—lpk 
j.k ^con^jk ^ohs /gx 

where Pcon(^) is the convolved theoretical power spec- 
trum, and -Pobs(^) is the observed power spectrum with 
the covariance matrix C. 



2. Model B: Halofit 

The Halofit model was developed by Smith et al. [15j . 
calibrated from simulations made by the Virgo consor- 
tium. It is an improvement on the then commonly used 
Peacock-Dodds formula , as it allowed for the modifi- 
cation of dark matter halos by continuing mergers. It is 
still only a simple empirical scaling however, and doesn't 
attempt to include the effects of redshift-space distortions 
or galaxy bias. 

We use the CAMB Halofit module to calculate the non- 
linear matter power spectrum Phf(fc) at redshift z. The 
bias is treated the same way as for the linear model. 



3. Model C: Simulation inspired fitting formulae 

Modelling the non-linear redshift space distortions of 
the galaxy power spectra is difficult, and one approach 
is to use large scale TV-body simulations. In Jennings 
et al. [H], the authors studied such effects in their sim- 
ulations, and produced fitting formula for the density, 
velocity divergence and cross-power spectra in the non- 
linear regime that provides a good fit to the simulations 
(up to fc ^ 0.2/iMpc~^). We briefly outline our interpre- 
tation/application of these formulae as follows. 

The angle averaged redshift-space power spectrum of 
matter can be written as |30j 

P„(fc,z) = A^{k)b^{z)Pss{k,z) (7) 
+2A2{k)fbPse{k, z) + A4k)f^Pge{k, z) , 

where Pgsik) is the non-linear density power spectrum, 
Peeik) is the power spectrum of the velocity divergence 
field often referred to as the velocity power spectrum, 
and Pse{k) is the cross-power spectrum. We use Halofit 
to calculate Pss{k, z). The velocity and cross-power spec- 
trum are calculated using the fitting formulae of Jennings 
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et al. |12i Eq. (15)] 4- Model D: Simulation inspired fitting formula with zero 

pairwise velocity 

P (k z = 0) = ^oVPssjk.z^ 0) + aiP^gjk, z = 0) 

' ' a2 + a^Pssik, z = 0) ' The fitting formulae used in model C were derived for 

(8) dark matter particles and not for halos. Setting all galaxy 
where Pss{k, z = 0) is the non-linear matter power spec- velocity dispersions to zero in model C provides a better 



trum at z = 0. The values of a are given in Table [III] for fit to the GiggleZ halo catalogues, and we have treated 
(xy) — (09, 59). this special case as a separate model. 



TABLE III: Best fit values of a in Eq. [Sjfrom Jennings et al. 



5. Model E: Empirical redshift space distortion damping 





Pee 


Pes 




-12462.1 


-12288.7 


ai 


0.839 


1.43 


02 


1446.6 


1367.7 


03 


0.806 


1.54 



The Pee{k) and Pse{k) must be scaled to redshift [T^ 
Eq. (17)], 

^^^-^^^)- D{z)+D^{z) + D^{z) ' 

where D is the growth factor. For a flat universe with a 
cosmological constant, the growth factor can be approx- 
imated as 



D{z) 



H{z) r°°dz'{l + z') 



Hn 



m{z') 



dz'{l + z') 
m{z') 



(10) 



where H{z) is given by Eq. [i] The power spectrum at 
redshift z is then given by [12, Eq. (17)] 



P,y{k,z = Q)-Pss{k,z = Q) 



^xy {k^ z) 

The coefficients. An, are calculated as 



+ Pssik,z) 
(11) 



An{k) = 



(12) 



where = fc-z is the cosine of the angle between the wave 
vector, k, and the direction of the line of sight, z. The 
growth rate, /, is given by / = dluD/dlna « n„^{z)'^'^^ 
with a being the scale factor, and the one-dimensional 
velocity dispersion, a^, given by 



2^1 f <Pk' Pee{k') 
~ 3 i (27r)3 



(13) 



Pgg is assumed to be spherically symmetric so the integral 
becomes 



2 1 



dk'Peeik'). 



" 3(27r)2 
We marginalise numerically over b 



(14) 



Non-linear structure formation will lead to increased 
peculiar galaxy velocities at low redshift (sometimes de- 
noted non- linear velocities), leading to a damping in the 
observed structure, which can be described by the em- 
pirical model [29] : 



Pgal(fc) = blP„,{k) 



Pm {k) 



-26./ 



1 + (kayfi)^ 
dfi 



d^ 



1 + {ka^^iY 



(15) 
(16) 



H- (kayfj.) 



1 + {ka^^iY 



where /i is the angle to the line of sight (defined as in 
Eq. 12), and tTi, is the pairwise galaxy velocity dis- 
persion in units of Mpc/ft,. Here the damping due to 
the pairwise velocity dispersion is a Lorentzian, rather 
than a Gaussian. The last equality is obtained using 
/3 = f /hr with / = Vlm{z)^'^^ ■ The velocity dispersion 
can then be calculated theoretically (e.g., from Eq 
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with Pgg{k) — f'^Pss{k))- The standard in the literature, 
which we follow here, is to set Pge{k) — pP\in{k)- 

Setting ay = OMpc//i we recover the linear bias case 
given in Eq. ]5] with the linear bias related to by by a 

simple numerical factor, b = 6. ^/^^l +7,^)^d,. Here 
we numerically marginalise over b^- 



6. Model F: N-body simulation calibrated 

All the non-linear effects are present in the GiggleZ 
simulation for a fiducial cosmology. Apart from the non- 
linear structure formation which we get from Halofit, we 
assume that any other effect is either negligible or can 
be scaled with cosmology. This is a simplified version of 
the method presented in Reid et al. P5], who considered 
three modifications of the linear power spectrum, treated 
independently. These modifications are: the damping of 
the BAO, the change in the broad shape of the power 
spectrum because of non-linear structure formation, and 
the bias because we observe galaxies in haloes in red- 
shift space rather than the real space matter distribution. 
Here we simplify it by substituting the complex scale- 
dependent bias of the Luminous Red Galaxies (LRGs) 
with the linear bias of WiggleZ galaxies. 
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For the BAO damped power spectrum we use the fol- 
lowing form 



^damped I 



ptrial / 
Mill 



(fc)/damp(fc)+/';,r(^)(l 



/damp) 1 

(17) 

where Pun is the linear power spectrum, and P^w is a 
'no-wiggles' power spectrum. Here the damping f act or is 
/damp(fc) = exp(-(fccri,)2), with ay given by Eq. 
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and 

Pee{k) = P\\n{k). The no- wiggles power spectrum (with- 
out the acoustic peaks) we get by spline-fitting the output 
power spectrum from CAMB following the approach of 
Reid et al. [55] and Swanson et al. [5^. The exact posi- 
tions of the nodes do not affect the final result as long as 
they are consistent for all sets of trial parameters. 

For the change in broadband power due to non-linear 
structure growth we use the Halofit model, applied as a 
ratio to the no- wiggles power spectrum (since Eq. 17 has 
already taken into account the effect on the BAO peak of 
non- linear structure formation). This ratio term is given 
by 



fhf 



■Phf,nw(fc) 



(18) 



Finally, we need to model the bias that comes from 
observing galaxies in halos. For this we make use of the 
GiggleZ simulations. -Pcfsgiozl^) given by a fifth-order 
polynomial fit to the power spectrum of a set of halos 
in the GiggleZ simulation chosen to match the clustering 
amplitude of the WiggleZ galaxies (details of this fitting 
formula can be found in Appendix lAl) , 



Dfid 
GiggleZ 



(k) = Ppoly(fid)(fc)Prmpcd(fc)/^„w(fc) • (19) 

where "fid" refers to the fiducial cosmology for the Gig- 
gleZ simulation. 

Combining all these terms together, we get (for a given 
trial cosmology) 



ptrial 
^gal 



(k) 



^2ptHal ,.^^M;.t(fc) Pl^(fc) ^afgglez(fc) 



damped 



(fc) 



ptrial 
nw 



and re-aranging we get 



(k) P^f^n„{k) Pdampcd(^) 

(20) 



^gTr'(fc) = b'PTn^ik) 



trial 



ptrial pfid 
^damped ^polyV'^7 

P^i^ik) Pfi'!nw(fc) 



(21) 



where is the linear scaling (related to galaxy bias) . We 
marginalize analytically over b. The second factor in Eq. 
[21] represents the smooth power spectrum of the trial 
cosmology. The third term defines the acoustic peaks 
and their broadening caused by the bulk flow motion of 
galaxies from their initial positions in the density field, 
and the fourth factor describes all the additional non- 
linear effects from the A'^-body simulation. 



the simulation and modelling that my not be present in 
WiggleZ, or we have included and that were not included 
in ? , Reid:2010] One particular case we tested was to 
see if introducing a strong damping may in fact lead to 
biased parameter values. Removing all damped power 
spectra in model F (so (T„ = and /damp(fc) = 1 in Eq. 
171 simplifies it to: 



p*rr^(fc) = 6^p*r'(fc) 



trial / 



^poly(fc) 

P^fik) 



(22) 



We marginalize analytically over b. 

We found this model provides a better fit to the Gig- 
gleZ halo catalogues, recovering the input fiducial cos- 
mology more accurately than in model F, and we have 
treated this special case as a separate model. 



B. Bias 

The GiggleZ simulations show that the effect of scale 
dependent bias is less than 1% for haloes with large-scale 
clustering amplitude matching that of the WiggleZ sam- 
ple, over the considered fc-interval, and consequently we 
only consider scale-independent bias [35] . Where possible 
we analytically marginalise over this bias, and otherwise 
numerically, as denoted for each of the approaches. 

The bias is required to be identical for all windows at 
the same redshift, so when we marginalise we assume a 
single parameter that specifies the bias at that redshift. 
We allow the bias to vary between redshifts, since the 
survey colour and magnitude cuts cause the galaxy pop- 
ulations observed in our sample to be different at different 
redshifts. 



C. A comparison of the approaches 

The power spectra for a fixed cosmology for models 
A-G are shown in Fig. [3| The different models tend to 
diverge for large k ik > 0.2 ft,Mpc~^), and these different 
predictions will have a significant effect on fitting the 
cosmological parameters. 

We tested the different approaches by fitting to the z — 
0.6 power spectrum of the GiggleZ simulation'^. For this 
we used two sets of parameter grids: versus (where 
/{, — r^b/f^m is the baryon fraction), and ^Im versus 
both with the remaining parameters fixed to the GiggleZ 
fiducial cosmology values. We chose these grids because 
the parameters are tightly related to the shape of the 
power spectrum. In both cases we obtain very similar 
conclusions so here we only present the results of 
versus /;,. 



7. Model G: N-hody simulation calibrated without damping 



The approach in model F was derived for a particular 
data set [28] . with several assumptions that went into 



^ This comparison was summarised in Riemcr-S0rensen et al. 1111 . 
but we repeat it here for the purposes of ongoing discussion, and 
the inclusion of model F. 
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FIG. 3: (color online), {top): A weighted average of the Wig- 
gleZ power spectra in the survey regions and redshifts, and the 
seven models described in section llV Al for the best fit cosmol- 
ogy of the best model (G). [bottom): The same models, but 
now the power spectrum has been rated to the linear predic- 
tion. In both plots the shaded region from k = 0.02 h}Apc~^ 
to k = 0.3/iMpc~^ is the range of data we are using in our 
analysis. This shows the divergence of models at large k and 
the reason why careful modelling is necessary. 



The ability to recover the input model is shown in Fig. 
|4| For fcmax < 0.2 /iMpc^^ most of the models produce a 
good fit, whereas for fcmax = 0.3 /iMpc~^ model B, C and 
E give reduced values above 1.5. Models B & C were 
calibrated to be good fits up to A; = 0.2/iMpc~^, so this 
is not surprising. The upper panel of Fig. [s] shows the 
for the fiducial GiggleZ cosmological parameters, which 
is a measure of how well the models recover the input pa- 
rameters. The lower panel shows the difference between 
X^ of the GiggleZ values and the best fit, indicating how 
far the best fit is from the input values. We assume that 
the A^-body simulation, which provides a complete census 
of the relevant non-linear effects, yields the most accurate 
clustering model. In this sense the good performance of 
model G (figures |4] and [s]) is a consistency check, and the 
variations of results produced by the other models are 
due to the breakdown in their performances compared to 
the simulation. 

It is clear that model G demonstrates the best "fit and 
recover" performance, and consequently we have chosen 
this model to be the default model of the WiggleZ Cos- 
moMC module. 
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FIG. 4: (color online). 1 and 2a contours for in the {Qm, 
fb} plane for fits to the GiggleZ power spectra using model 
A-G. All other cosmological parameters are held fixed at the 
GiggleZ fiducial cosmology (WMAP 5-year best fit values). 
The GiggleZ parameter values are marked by a cross. The red 
ellipses (model E) lie directly on top of the green ones (model 
B), and so are not visible. This shows that the choice of 
model has a significant effect on the cosmological constraints 
from the power spectrum. Our preferred model is the Model 
G, which successfully reproduces the input parameters of the 
GiggleZ simulation. 
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FIG. 5: (Top) Reduced of models fitted to the A-body 
simulation halo catalogue for the GiggleZ fiducial cosmology 
values. In absence of systematic errors the models should 
recover the input cosmology with x^/dof ~ 1. (Bottom) Dif- 
ference in reduced x^ values when using the GiggleZ fiducial 
cosmological parameters and the best fit values. 



D. Systematics 

In measuring the power spectrum from the WiggleZ 
data a number of choices were made in modelling cer- 
tain unknown factors, which may impact on the result- 
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ing data product and cosmological parameters. These 
were the nature of the radial selection function, N{z), 
and the fiducial cosmological parameters used to gener- 
ate the power spectra. We test a number of possible sys- 
tematics to see how much they impact on the resulting 
cosmological parameter constraints. 

Firstly there is the choice of radial selection function, 
N{z). The number of galaxies in the WigglcZ sample 
will vary not just as a function of angular position on the 
sky, but also as a function of redshift, because the sur- 
vey colour and magnitude cuts cause the observed galaxy 
populations to be different at different redshifts. It is dif- 
ficult to calibrate this directly, so we assume some smooth 
function form for N{z). We fit the same N{z) function to 
each region of the sky, but a different one for each observ- 
ing priority band (WiggleZ observations were prioritised 
by magnitude, with apparently-faint galaxies given high- 
est priority and each band representing an equal interval 
of the range 20.0 < r < 22.5). The default choice for 
the function is a Chebyshev polynomial, but alternatives 
include a double-peaked Gaussian, and a cubic-spline fit. 
In Fig. |6]we show the effect of these choices on the recov- 
ered joint cosmological constraints on the physical matter 
density ilmh^ and the baryon fraction (keeping the re- 
maining parameters fixed at the WMAP7 best fit values). 
We see the effect on the recovered cosmological param- 
eters from the choice of N{z) is small, with the error 
ellipses being consistent between each other at the two- 
sigma level. This is as expected given that the data are 
identical. Since there is no obvious large systematic, our 
choice of N(z) is motivated by our initial confidence in 
which of these models will work best. The double Gaus- 
sian has too few parameters to describe real variations, 
whereas the spline fit has too many parameters and may 
over fit, removing real fluctuations. The Chebyshev poly- 
nomial (whose order is truncated as driven by the data) 
is a compromise between these effects. 

Secondly there is the choice of fiducial cosmological pa- 
rameters in converting the redshifts to distances. We in- 
clude the Alcock-Paczynski scaling factor when comput- 
ing the predicted power spectrum in the likelihood calcu- 



lation, as given in section III We check the effect of as- 
suming the wrong cosmology by re-measuring the power 
spectra using a fiducial matter density of flf^ = 0.37. 

From Fig. [6] it is clear that the choice of fiducial cos- 
mological parameter values has a small effect on the best 
fit values, though the size of the error ellipse is expanded. 

There is also the choice of the range of fc-values to 
fit the data to. Throughout our discussion of these 
models, and further parameter analysis, we have fixed 
the lower limit of power spectrum fitting to be k^in = 
0.02 h Mpc""'^, which corresponds to the largest modes ob- 
served in each of the WiggleZ regions (the final results are 
not very sensitive to the exact value). The upper limit 
on the fitting range was fixed to kmax = 0.3/iMpc^^ 
after running a series of flat ACDM fits to WMAP -f 
WiggleZ varying fcmax- As seen in Fig. [7j the uncer- 
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FIG. 6: (color online). 1 and 2a contours for in the 
/{,} plane for fits to the WiggleZ power spectra, checking for 
systematic errors using a number of different data analyses. 
The red dotted lines (spline fit) lie directly on top of the black 
(Chebyshev), and so are not visible. 



start to increase again (due to the effect of system- 
atics), before decreasing to roughly the same value at 
fcmax = 0.5 /i Mpc~^. Furthermore, the mean values of 
the parameter probability distributions are roughly con- 
sistent in the range of maximum fc considered, similar to 



what we saw in section IV C when making comparisons 
to simulations. Our choice of fcmax = 0.3 h Mpc~^ is then 
a compromise between minimising the parameter errors, 
and trusting that our simulations are accurate down to 
those small scales. 




FIG. 7: (color online). Measurement of flmh^ (top) and /b 
(bottom) as a function of fcmax for the WiggleZ data, with 
all other cosmological parameters fixed to the WMAP 7-year 
best fit values, with those two being measured simultaneously. 



tainties decrease until k„ 



0.3/iMpc^^ and then 



11 



E. WiggleZ power spectrum fit 

To check the consistency between the cosmological pa- 
rameters derived from WMAP and WiggleZ, we fitted the 
WiggleZ power spectra alone varying the physical matter 
density, ftmh^, and the baryon fraction, /;,, keeping the 
remaining parameters fixed at the WMAP7 best fit val- 
ues. The results are shown in Fig. [8] where we also show 
the constraints from the CMB (WMAP 7-year) in com- 
parison. The two data sets give consistent cosmological 
parameters in this model. There is a tendency for Wig- 
gleZ alone to favour slightly smaller values of the baryon 
fraction in comparison to the CMB data alone, but not 
at any statistically significant level. 



Cosmological Parameter 



symbol 
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FIG. 8: (color online). The probability contours in the 
{Q,mh? , fb} plane for WiggleZ alone (black solid line), with 
all other cosmological parameter fixed to the WMAP 7-year 
best fit values. We also plot the contours from the WMAP 
data for a ACDM model, marginalising over all other parame- 
ters (blue filled contour). We see good agreement between the 
two data sets, though the WiggleZ data does favour slightly 
smaller values of the baryon fraction. 



V. COSMOLOGICAL CONSTRAINTS 

In this section we will discuss the different cosmological 
parameters and combinations of parameters (models) to 
be tested, and present the results of our analysis. 



A. Cosmological Models 

The approach we take in this paper is to start with the 
flat ACDM model as our benchmark. We first check the 
consistency of our new WiggleZ P{k) data with previous 
data sets, within the flat ACDM model. The basic model 
has flve parameters, as described in the introduction and 
summarised in table HVl 



Base parameters (for Flat ACDM) 
Physical baryon density 
Physical CDM density 
Angular size of the sound horizon 
at decoupling 

Optical depth of reionization 
Scalar spectral index 
Amplitude of the scalar perturbations 
Sunyaev-Zel'dovich (SZ) template normalization Asz 



T 

Us 

A., 



Extension parameters 

Dark energy Equation of state w 

Curvature f2fc 

Neutrino density fraction fu 

running of the spectral index rirun 

Tensor to scalar amplitude ratio r 
Derived parameters 

Matter density fi™. 

Dark energy density JIa or SIde 

Hubble parameter Ho or h 

Matter power dispersion at 8/i~^Mpc erg 

Combined mass of the neutrino species Yl 

TABLE IV; Cosmological parameters to be varied. 



Section Model 



Extra Parameter (s) 



VDl 


Flat ACDM 


none 


VD2 


Flat wCDM 


w 


VD3 


ACDM 




VD4 


wCDM 




VD5 


massive neutrinos 


.u 


VD6 


running spectrum 


^run 


VD7 


gravitational waves 


r 



TABLE V: Cosmological models to be analysed. 



To extend the ACDM model, we consider parame- 
ters whose existence is required or suggested by either 
external datasets (for example neutrino oscillation ex- 
periments demonstrating that the neutrinos have mass) 
or from theoretical predictions (e.g., primordial tensor 
perturbations generated by inflation). While there is a 
wealth of possible extension parameters to choose from, 
we limit ourselves to these five: the equation of state of 
the dark energy w (assumed to be a constant) , the curva- 
ture, the running of the primordial density power spec- 
trum, the tensor to scalar ratio, and the density fraction 
of massive neutrinos. 

The standard set and extension parameters are given 
in table IIVI The different models we will consider will 
consist of the base parameter set plus one or more extra 
parameters, by themselves or in combination with others. 
These are shown in table |Vl 

If a parameter is not varied in some model, and can- 
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not be derived from the other parameters that are being 
varied, then it is held fixed at some defauh value. These 
default values are zero, except in the case of the equa- 
tion of state of the dark energy, where the default value 
is m; = — 1 (i.e. a cosmological constant). 



B. Parameters and the physics being constrained 

The cosmological observations that we make constrain 
different physics, and so it is important, especially with 
many-parameter models, to understand exactly which 
data are contributing to constraints on which parame- 
ter or parameters. The major cosmological datasets are: 

CMB: The CMB is probably the most powerful tool 
available to a cosmologist. The features of the CMB 
power spectrum are the Sachs- Wolfe plateau which mea- 
sures As, and gives weak constraints on the dark energy 
parameters floE and w through the Integrated Sachs 
Wolfe effectl the position of the acoustic peaks which 
function as a standard ruler measuring flm, ^^A, and Hq] 
the relative height of the acoustic peaks which measure 
flm and fib! and the damping tail, which gives some weak 
constraint on r. There is also the polarisation power 
spectrum, which gives stronger constraints on r. Finally, 
as the CMB is measuring the amplitude of density fluc- 
tuations over a range of scales it provides a constraint on 
the spectral index Ug ■ For the analysis presented here we 
have used the CMB observations from the WMAP 7-year 
data release [3T] . 

BA 0: The Baryon Acoustic Oscillations are a standard 
ruler, and so give constraints on those parameters that 
relate to distance, fi^, r^DE, and w. Their size is set in 
co-moving coordinates, so it is not possible to use them 
to measure the Hubble parameter today (Hq) without 
reference to some physical size, such as from the CMB. 
Despite their name, the size of the ruler depends on the 
baryon density only weakly, and since our measurement 
of ilhh^ from the CMB is already very accurate, this 
parameter is usually fixed when using BAO data. For 
the analysis presented here we have used the BAO scale 
as measured by SDSS [23132]. 

Note that we do not include the WiggleZ BAO mea- 
surements, as they come from the same galaxy sample 
as our P{k) measurement, and so we are being conser- 
vative in the light of potential covariance. We compare 
the BAO results from WiggleZ with the results from the 
WiggleZ P{k) in a couple of cases of interest (curvature 
and the dark energy equation of state), and find fairly 
good agreement in the cosmological constraints from the 
two datasets. 

LSS: The large-scale structure of matter gives informa- 
tion about the matter density flm through the broadband 
shape of the power spectrum, with some modulation by 
baryons flh or neutrinos and the amplitude of fluc- 
tuations as (though this measurement is somewhat de- 
generate with the linear bias). It also contains the BAO 
peaks, and so in principle can be used to give the same 



distance information as BAO. However, when measuring 
the full power spectrum the overall shape is more im- 
portant than the wiggles, so the BAO features have only 
a small effect on the cosmological constants for the full 
power spectrum. The LSS data used here is SDSS LRG 
sample [25], and the WiggleZ P{k) measurement in this 
paper. (Note that again, because of previously mentioned 
covariance between the two data sets, we never combined 
the SDSS LSS and BAO data in the same analysis). In 
Fig. [9) we show the effect of varying the matter density, 
baryon fraction, spectral index and neutrino fraction on 
the matter power spectrum. 

SNIa: Type-la Supernovae are standard candles, and 
so also give constraints on the parameters that control 
the cosmological distances, ^^de and w. Their ab- 
solute magnitude is not known, and is degenerate with 
the Hubble parameter, so this combination is normally 
marginalised over when using this data set. We have used 
the SNLS data set 33-35J. 

Hq: Direct measurements of the Hubble parameter 
from local standard candles or rulers also play a key role, 
establishing a normalising 'ruler' in the cosmological dis- 
tance ladder. Here we have used the Hq data point given 
in l36l. 



C. Analysis method 

The most commonly used approach in cosmology is 
a Bayesian analysis of the the cosmological data. By 
Bayes' theorem, the posterior probability distribution of 
a parameter, 9, given a model, M, and a dataset, D, is 
given by 

^mM) = ^^^MS^= (23) 



P{D\M) 



E 



where P{9\M) (or 7r(6')) is the 'prior' parameter distribu- 
tion, and P{D\e,M) (or C) is the 'likelihood*.' The dif- 
ferent models we test are given in Table |Vj In this paper 
we focus on the posterior distributions (parameter con- 
straints) of the different models for different data combi- 
nations, and set aside comparisons between the models, 
except in view of their consistency with the data. 

For the parameter priors, we assume a uniform prob- 
ability distribution, with limits given in Table |VI[ These 
prior ranges are chosen such that they are much larger 
than the expected posterior ranges (with certain excep- 
tions, e.g., neutrino mass), so that we do not not miss 
regions of interest in the parameter space, and do not 
necessarily correspond to any physical intuition. 



* The normalising coefficient P{D\M) (or E) is the 'model likeli- 
hood' or 'evidence', and can be used to evaluate the probability 
of the model as a whole (for examples of this use in cosmology, 
see |37P3) V 




FIG. 9; (color online). The matter power spectrum plotted for different values of the cosmological parameters. For each plot 
we fix the other parameters to the fiducial cosmological values while varying only that particular one. The colour bar indicates 
the value of the parameter. The dashed black curve is the power spectrum for the fiducial cosmology, and is the same in all 
four plots. 



Due to the complexity of the datasets under considera- 
tion it is impossible to solve for the individual parameter 
posteriors analytically. Instead we must proceed by sam- 
pHng of the parameter space, and averaging ('marginali- 
sation') of the likelihoods over the parameter ranges. For 
this we used the CosmoMC cosmological analysis package 
[TT] . CosmoMC is a Markov Chain Monte Carlo pack- 
age that compute the likelihoods of all the cosmological 
datasets we use. We sample the parameter space using 
the Metropolis-Hastings algorithm, which comes as stan- 
dard in the package. 



We have created an additional WiggleZ power spec- 
trum module for the code, that includes the modelling 
details of Model G. The module computes the likelihood 
of each WiggleZ redshift bin independently, but computes 
contributions from all regions simultaneously. It does this 
in order to marginalise over the linear bias b, using Eq. 
[6j assuming the bias is the same for every region at a 
given redshift. We make the power spectrum data and 
the module available with this paper. 



D. Results 



FlatACDM 



We initially investigate how combining the WiggleZ 
power spectrum data with the other cosmological data 
sets constrains parameters in the concordance ACDM 
model. The WiggleZ power spectrum data alone are 
not able to constrain all the parameters, so we add the 
WMAP 7-year CMB data. We then considered the effect 
of adding extra datasets to the CMB -I- WiggleZ pairing. 

The marginalised constraints are given in table |VII| 
The two-dimensional parameter constraint contours for 
{Qrrn Hq} are shown in Fig. 10 and {rts, cts} are shown 
in Fig. [nj 

We find that adding the WMAP CMB data to the 
WiggleZ P{k) data considerably shrinks the size of the 
confidence contours, as well as allowing us to constrain 
other parameters we would not be able to using only the 
WiggleZ data (such as r). The contours are smaller than 
both individual datasets by themselves, and the Wig- 
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Parameter 


minimum 


maximum 






0.005 


0.34 


0.90 




0.01 


1.5 




10061 


0.5 


10 


0.85 




0.01 


0.5 




ris 


0.5 


1.5 




log(10^''A,) 


2.7 


4 


0.80 


Asz 





2 




ID 


-3 













0.75 


-0.3 


0.3 







1 




^run 


-0.2 


0.2 


0.70 


r 





1 





TABLE VI: Prior ranges on the cosmological parameters. 
We assume uniform probability distribution between these 
ranges. 
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FIG. 10: (color online). The two-dimensional constraints in 
the {fim, Ho} plane for a flat ACDM model, using a num- 
ber of different datasets (see legend). The WiggleZ data are 
consistent with the CMB data (from WMAP), reducing the 
size of the error ellipse, but pushing the preferred value of Q.m 
upwards. 
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FIG. 11: (color online). The two-dimensional constraints in 
the {(Tg, nji} plane for a flat ACDM model, using a number of 
different data combinations (see legend). We find that com- 
bination of LSS data (from either WiggleZ or SDSS) with the 
CMB data (WMAP) only improves the constraints on the 
spectral index ris marginally, but improves the constraints on 
the amplitude of fluctuations ag by a large amount. 



2. FlatwCDM 

Here we consider models where the equation of state 
of the dark energy, w, is allowed to be different from — 1 
but kept constant with time. Again we only consider flat 
models. 

In addition to CMB -I- WiggleZ, we also compare con- 
straints on the constant equation of state w using the 
WiggleZ P{k) with the SNLS Type-la Supernova compi- 
lation, and then combining the two. The supernova data 
are far more powerful at measuring distances, having 
many more data-points spread over a range of redshifts, 
and so perform much better than the WiggleZ P(fc), as 
shown in Fig. 
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The marginalised constraints are given in table |VIll 



3. ACDM 



gleZ data are complimentary to WMAP, narrowing the 
confidence contours of the total matter density around 
rim = 0.29 (and since the distance to recombination 
is well constrained in this model, preferring values of 
h = 0.689). In terms of the parameters that govern the 
primordial power spectrum, n^, and a^, we see from Fig. 
11 that the WiggleZ data improves constraints only on 
the amplitude of fluctuations, but does not help with the 
spectral index (a similar effect is seen using CMB -I- LSS 
data). This is due to the fact that the constraint from 
WMAP on Us is already very good, and the large scale 
measurements of the matter power spectra at late times 
do not offer enough accuracy over a large enough range 
of scales to be effective. 



Here we consider models where the curvature is allowed 
to be different from flat, allowing fife to be non-zero. The 
marginalised constraints are given in table [VTll The two- 
dimensional parameter constraint contours for {flm, ^a} 
are shown in Fig. [13] 

We find two, slightly separated, high-likelihood peaks. 
One of these is close to flatness, with a matter density 
ilm = 0.29 and Hubble parameter h = 0.69. The other 
is slightly further away, with matter density VL^ = 0.52 
and Hubble parameter h — 0.50. These are allowed by 
the CMB due to the extra dynamical freedom obtained 
by letting the Universe move away from flatness. The 
high matter density peak is not ruled out by including 
the WiggleZ data, as both of these parameter sets have 
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Model 


Parameter 




CMB + WiggleZ 


± Ho 


+ SN-Ia 


+ BAG 


+ J?o ± BAG 


Flat ACDM 






2.238 ±0.052 


2.255 ±0.050 


2.240 ±0.053 


2.239 ±0.050 


2.253 ±0.050 








0.1153 ± 0.0027 


0.1145 ±0.0026 


0.1150 ±0.0028 


0.1152 ±0.0024 


0.1146 ±0.0024 




woe 




1.039 ±0.002 


1.040 ± 0.002 


1.039 ±0.003 


1.039 ±0.002 


1.039 ± 0.002 




T 




0.083 ±0.014 


0.084 ±0.014 


0.083 ±0.014 


0.083 ±0.014 


0.084 ±0.014 




Us 




0.964 ±0.012 


0.968 ±0.012 


0.965 ±0.013 


0.964 ±0.012 


0.968 ±0.011 




log(10^°yl. 


) 


3.084 ±0.029 


3.086 ± 0.029 


3.085 ± 0.030 


3.083 ±0.029 


3.086 ± 0.029 








0.290 ±0.016 


0.283 ±0.014 


0.288 ±0.017 


0.289 ±0.013 


0.284 ±0.012 




/fo[kms-i 


Mpc-i] 


68.9 ± 1.4 


69.6 ± 1.3 


69.1 ± 1.6 


69.0 ± 1.2 


69.5 ± 1.2 








0.825 ±0.017 


0.825 ±0.017 


0.825 ±0.017 


0.825 ±0.017 


0.825 ±0.017 


Flat wCDM 






2.265 ±0.062 


2.253 ±0.057 


2.228 ±0.055 


2.247 ±0.056 


2.253 ±0.056 








0.1164 ±0.0036 


0.1146 ±0.0030 


0.1157 ±0.0030 


0.1147 ±0.0029 


0.1148 ±0.0030 




ime 




1.039 ± 0.003 


1.039 ±0.003 


1.038 ±0.003 


1.039 ±0.003 


1.039 ± 0.003 




T 




0.084 ±0.015 


0.084 ±0.014 


0.082 ±0.014 


0.084 ±0.014 


0.084 ±0.014 








0.975 ±0.019 


0.968 ±0.014 


0.962 ±0.014 


0.967 ±0.014 


0.968 ± 0.014 








3.096 ±0.031 


3.086 ± 0.030 


3.082 ±0.029 


3.085 ± 0.030 


3.086 ± 0.030 




w 




-0.525 ±0.293 


-1.007 ±0.084 


-1.062 ±0.072 


-0.973 ±0.086 


-1.008 ±0.085 








0.487 ±0.132 


0.283 ±0.018 


0.844 ±0.028 


0.294 ±0.018 


0.284 ±0.018 




Ho 




55.2 ±8.4 


69.7 ± 2.1 


70.5 ±2.3 


68.4 ±2.0 


69.7 ±2.1 








0.664 ±0.110 


0.826 ± 0.032 


0.844 ±0.028 


0.815 ±0.033 


0.827 ±0.032 


ACDM 


100f7b/i^ 




2.215 ±0.055 


2.263 ± 0.054 


2.256 ±0.054 


2.252 ±0.054 


2.262 ±0.052 








0.1118 ±0.0039 


0.1162 ±0.0039 


0.114 ±0.0042 


0.1150 ±0.0038 


0.1161 ±0.0038 




1000 




1.038 ±0.003 


1.040 ±0.003 


1.040 ±0.003 


1.040 ±0.003 


1.040 ±0.003 




T 




0.086 ±0.014 


0.088 ±0.015 


0.089 ±0.014 


0.088 ±0.015 


0.088 ±0.014 




Us 




0.958 ±0.013 


0.970 ±0.013 


0.969 ±0.013 


0.968 ±0.013 


0.969 ±0.013 




log(10'''yl. 


) 


3.072 ±0.031 


3.101 ±0.031 


3.096 ±0.031 


3.096 ±0.031 


3.101 ±0.030 








0.454 ±0.058 


0.287 ±0.029 


0.303 ±0.038 


0.302 ±0.020 


0.288 ±0.016 








-0.046 ±0.017 


0.001 ±0.008 


-0.005 ±0.012 


-0.004 ±0.006 


0.000 ±0.005 




Ho [kms-i 


Mpc-i] 


54.65 ±3.8 


69.86 ±3.6 


67.7 ±4.7 


67.6 ±2.3 


69.9 ±3.6 




0-8 




0.782 ±0.024 


0.838 ± 0.023 


0.825 ±0.026 


0.829 ±0.022 


0.838 ± 0.023 


wCDM 


100f7b/i^ 




2.231 ± 0.058 


2.248 ±0.059 


2.205 ±0.055 


2.244 ±0.057 


2.246 ± 0.056 








0.1117 ±0.0038 


0.1133 ±0.0045 


0.1094 ±0.0041 


0.1130 ±0.0041 


0.1131 ±0.0040 




looe 




1.039 ± 0.003 


1.039 ± 0.003 


1.038 ±0.003 


1.039 ±0.003 


1.039 ± 0.003 




T 




0.086 ±0.015 


0.084 ±0.014 


0.083 ±0.014 


0.084 ±0.014 


0.084 ±0.014 








0.964 ±0.015 


0.966 ±0.014 


0.956 ±0.013 


0.966 ±0.014 


0.966 ±0.014 




log(10i"^. 


) 


3.075 ±0.032 


3.080 ± 0.034 


3.054 ±0.032 


3.079 ± 0.032 


3.079 ±0.031 




w 




-0.502 ±0.160 


-1.102 ±0.219 


-1.215 ±0.117 


-1.060 ±0.177 


-1.113 ±0.169 








0.614 ±0.061 


0.283 ±0.018 


0.354 ±0.041 


0.289 ±0.020 


0.278 ±0.017 




fife 




-0.052 ±0.018 


-0.005 ±0.012 


-0.032 ±0.015 


-0.005 ±0.010 


-0.005 ±0.009 




Ho [kms-1 


Mpc-i] 


46.9 ±2.6 


69.4 ± 2.1 


61.3 ±4.2 


68.5 ±2.1 


69.9 ± 1.8 




erg 




0.654 ±0.051 


0.836 ± 0.040 


0.818 ±0.028 


0.826 ± 0.040 


0.842 ± 0.037 



TABLE VII: Cosmological parameter constraints combining WiggleZ data with other cosmological datasets, for different cos- 
mological models. 



the same broad-spectrum shape. We show this by plot- 
ting the relevant power spectrum for the two peaks in Fig. 
[141 Here we show both the initial linear power spectra 
and final prediction power spectra (convolved with the 
region-averaged window function at z = 0.6) in compar- 
ison with the data, and show that both models are good 
fits. However, if we were able to precisely determine the 
position of the turnover in the matter power spectrum. 



we would be able to rule out one of these models. (For 
further discussion of using the WiggleZ data to determine 
the turnover position and using the turnover to measure 
the matter density, 0^, see reference [M].) We can also 
rule out the high f2j„, fow Hq peak using other data, 
such as distance data from SN-Ia or BAO, or iJo mea- 
surements. This is shown in Fig. |13[ where the second 
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FIG. 12: (color online). The two-dimensional constraints in 
the {rim, plane for a flat wCDM model (see legend). The 
supernovae data (solid dark blue) is far more powerful at mea- 
suring distances, having many more data-points spread over 
a range of redshifts, but the addition of the WiggleZ data 
(thick red) has an effect on the best fit Q,ra value, pushing it 
slightly higher. 




FIG. 14: (color online). The matter power spectrum for two 
different sets of cosmological parameters, both of which are 
good fits to the ACDM model using CMB -I- WiggleZ data. 
The solid curves give the predicted linear P{k). The dashed 
curves give are the same P{k) after the non-linear effects and 
window functions have been applied, to compare to the data 
sit z — 0.6 (black error bars). 




FIG. 13: (color online). The two-dimensional constraints in 
the {^Im, ^a} plane for a ACDM model, for a number of 
different data compilations (see legend). The straight black 
line gives the location of Qm + = 1, i.e. a flat universe. 
The WiggleZ data are consistent with the CMB data (from 
WMAP) with a peak likelihood close to flatness, but have a 
second likelihood peak for larger values of Qm , and away from 
the flatness line. 



peak vanishes once the other data sets are applied. 

We also compare the constraints in the {0^, ^a} plane 
between the WiggleZ P{k) data and the WiggleZ BAO, 
as shown in Fig. [15] We see that only one of the high 
likelihood peaks in the P{k) result is favoured by the 
WiggleZ BAO measurement. 




FIG. 15: (color online). The two-dimensional constraints in 
the {flm, Qa} plane for a ACDM model, comparing the Wig- 
gleZ P{k) data with the WiggleZ BAO data (see legend). The 
straight black line gives the location of flm + ^a = 1, i.e., a 
fiat universe. The WiggleZ P{k) data are consistent with the 
BAO from WiggleZ, for the peak likelihood close to flatness, 
but not for the second peak at high f2,„ values. 



4. wCDM 

Here we consider models where both the equation of 
state of the dark energy, w, is allowed to be different from 
— 1 (but still constant with time) and also the curvature 
is allowed to be different from flat, allowing fi/j to be non- 
zero. Dark energy dynamics (a constant equation of state 
different from —1 or an equation of state which varies 
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FIG. 16: (color online). The two-dimensional constraints in the {Qk, w} plane (left) and {Qk, w} plane (right) for a uiCDM 
model, using a number of different datasets (see legend). The solid black cross gives the location of Flat ACDM, i.e. a flat 
universe with the equation of state w = — 1. 



with scale factor) has a somewhat degenerate signal with 
small curvature values 05]. If we were to detect w ^ 
— 1, it would be necessary to test both the equation of 
state and curvature simultaneously, to be sure that we 
had actually detected time-evolving dark energy (and so 
falsified the cosmological constant as the explanation for 
the late-time acceleration). 

The marginalised constraints are given in table |VII[ 
The two-dimensional parameter constraint contours for 
{rifc, w} and {rife. Ho} are shown in Fig. 16 We find 
that the CMB + WigglcZ data by themselves prefer a 
peak far away from the w = — 1, fife = values of flat 
ACDM. This peak is similar to the low Hq, high flm peak 
that appears in the analysis of the ACDM model, as it 
has the same similar best fit values for those parameters. 
As shown in Fig. 14 the WiggleZ P{k) data by itself 
cannot distinguish between these two sets of parameter 
values, and requires extra distance information in order 
to break the degeneracy. However, for this model, the 
ACDM point in the {Jlfe, w} plane is disfavoured at more 
than 2-sigma. While we could consider this as marginal 
detection of some new physics, it is more likely to be 
a coincidence, randomly occurring due to the number of 
different data sets and models under consideration in this 
paper. It could also be a product of some systematic er- 
ror, such as calibration of the selection function, popula- 
tion of haloes with WiggleZ galaxies in the GiggleZ sim- 
ulations, or non-linear model fitting. The option that it 
is a result of coincidence or systematic error rather seems 
especially true, given that most other datasets strongly 
disfavour this peak (the very low value of Hq is devoured 
by the Hq data for example). We discuss the goodness 
of fit improvement for this model, and consistency with 
other datasets, in section |VE] 

When the CMB -I- WiggleZ is combined with other dis- 
tance data, the constraints on both w and Ofe are only 
twice as large as when the parameters are considered indi- 



vidually. The best constraints on both parameters come 
from the CMB + WiggleZ -I- BAG -I- Hq combination, 
where w = -1.11 ± 0.17 and flk = -0.005 ± 0.009. This 
demonstrates that combining CMB, P{k) and other dis- 
tance data together provides reasonable constraints on 
these two extension parameters in combination. 



5. Massive neutrinos 

Recent particle physics experiments measuring neu- 
trino oscillations have show that mass differences exist 
between the different species [IHl HZ] ; but there is no cur- 
rent experiment that can measure the absolute neutrino 
mass. The massive neutrinos provide a form of 'Warm 
Dark Matter', suppressing structure formation on large 
scales. Measurements of large scale structure can con- 
strain the total neutrino density, and so the sum of the 
particle masses. Though we fail to directly detect the 
neutrino mass, we find the total mass of all three species, 
^m^, to be less than 0.58eV at 95% confidence when 
combining WiggleZ with the CMB. This upper limit de- 
creases to 0.34eV when we include distance data (BAG 
and Hq). 

Some of these results were first presented in Riemer- 
S0rensen et al. |llj using importance sampling of the 
WMAP MCMC chains. Here we present the output of 
the CosmoMC module and note that the two analyses 
give consistent results. The contours of the {Hp V nii,} 
and {^Im, planes are shown in Fig. |17[ The 

marginalised constraints are given in table [Villi 

It is worth noting how unconstrained rim becomes from 
WMAP -I- WiggleZ when allowing the neutrino mass to 
vary. This is because there are now effectively two types 
of Dark Matter (ordinary CDM and massive neutrinos 
acting as Warm Dark matter) , and the neutrinos can play 
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FIG. 17: (color online). (Left) The two-dimensional constraints in the {"^m,^, Q.m} (Left) and {J^^m^, Ho} (Right) planes for 
a massive neutrinos model, using a number of different data combinations (see legend). 



Model 


Parameter 




CMB + WiggieZ 


± Ho 


+ LSS 


+ BAG 


+ Ho + BAG 


ACDM + 


m„ lOOfifc/i^ 




2.240 ± 0.054 


2.264 ±0.051 


2.244 ±0.052 


2.250 ±0.051 


2.266 ±0.051 








0.1183 ±0.0052 


0.1157 ±0.0030 


0.1179 ±0.0031 


0.1161 ±0.0027 


0.1153 ±0.0025 




woe 




1.039 ±0.003 


1.040 ±0.002 


1.039 ±0.002 


1.039 ±0.002 


1.040 ±0.002 




T 




0.084 ±0.014 


0.086 ±0.014 


0.083 ±0.014 


0.085 ±0.014 


0.086 ±0.014 




Us 




0.964 ±0.013 


0.971 ±0.012 


0.965 ±0.012 


0.967 ±0.012 


0.971 ± 0.012 




/opflO^" 




3.089 ± 0.030 


3.090 ± 0.030 


3.088 ± 0.030 


3.087 ±0.030 


3.089 ±0.030 








0.042 


0.029 


0.034 


0.029 


0.026 








0.58 


0.32 


0.45 


0.38 


0.34 








0.316 ±0.042 


0.293 ±0.018 


0.311 ±0.021 


0.298 ±0.016 


0.291 ±0.014 




Ho [kms-i 


Mpc-1] 


66.9 ±2.7 


68.8 ± 1.5 


67.2 ± 1.7 


68.2 ± 1.3 


68.9 ± 1.2 








0.787 ±0.030 


0.798 ± 0.026 


0.791 ±0.027 


0.794 ±0.026 


0.799 ±0.025 


ACDM + 


rlrun WOQhh^ 




2.167 ±0.070 


2.237 ±0.064 


2.188 ±0.058 


2.206 ±0.055 


2.249 ±0.055 








0.1219 ±0.0056 


0.1185 ±0.0050 


0.1201 ±0.0035 


0.1181 ±0.0032 


0.1171 ±0.0035 




woe 




1.039 ±0.003 


1.040 ±0.003 


1.039 ±0.003 


1.039 ±0.003 


1.040 ± 0.002 




T 




0.090 ±0.016 


0.091 ±0.016 


0.090 ±0.016 


0.091 ±0.016 


0.092 ±0.016 




Us 




0.915 ±0.032 


0.951 ±0.029 


0.926 ±0.024 


0.935 ±0.022 


0.958 ±0.023 








-0.040 ±0.022 


-0.016 ± 0.021 


-0.033 ±0.018 


-0.027 ±0.017 


-0.011 ±0.018 




log(10^°A,) 


3.113 ±0.035 


3.113 ±0.035 


3.108 ±0.034 


3.106 ±0.034 


3.110 ±0.035 








0.332 ±0.039 


0.304 ± 0.030 


0.318 ±0.022 


0.305 ±0.018 


0.295 ±0.019 




Ho [kms-1 


Mpc"^] 


66.0 ±2.6 


68.2 ±2.3 


66.9 ± 1.7 


67.9 ± 1.5 


68.9 ± 1.6 








0.839 ±0.020 


0.842 ± 0.020 


0.835 ±0.018 


0.832 ±0.018 


0.840 ±0.019 


ACDM + 


r lOOfib/i^ 




2.268 ±0.063 


2.296 ±0.060 


2.257 ± 0.00056 


2.268 ±0.056 


2.288 ± 0.053 








0.1157 ±0.0036 


0.1142 ±0.0033 


0.1168 ±0.0029 


0.1158 ±0.0028 


0.1149 ±0.0027 




woe 




1.039 ±0.003 


1.040 ±0.003 


1.039 ±0.003 


1.040 ± 0.003 


1.04 ± 0.002 




T 




0.086 ±0.014 


0.089 ±0.015 


0.084 ±0.014 


0.086 ±0.014 


0.0880 ±0.014 




ris 




0.974 ±0.017 


0.982 ±0.016 


0.971 ±0.014 


0.974 ±0.014 


0.979 ±0.014 




log(10^°A. 


) 


3.098 ± 0.030 


3.101 ±0.030 


3.096 ± 0.030 


3.098 ± 0.030 


3.101 ±0.030 




r 




< 0.21 


< 0.23 


< 0.18 


< 0.20 


< 0.21 








0.289 ±0.023 


0.277 ±0.019 


0.297 ±0.018 


0.289 ±0.016 


0.282 ±0.015 




Ho [kms-i 


Mpc-i] 


69.3 ±3.6 


70.5 ± 1.9 


68.6 ± 1.61 


69.3 ± 1.5 


70.0 ± 1.4 




0-8 




0.835 ±0.018 


0.833 ±0.018 


0.837 ±0.017 


0.836 ±0.018 


0.835 ±0.018 



TABLE VIII: Cosmological parameter constraints combining WMAP ± WiggleZ data with other cosmological datasets, as- 
suming a flat ACDM model with some extra component. Uncertainties are la and upper limits are 95% confidence level. 
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some of the role in structure formation normally taken by 
CDM. This degeneracy is best broken by adding the BAO 
scale, as the extra distance data give better constraints 
on the combined matter density, and so reduce the size 
of the role neutrinos can play. 

The best constraints come from the CMB + WiggleZ 
+ Hq data compilation, and adding the BAO data on 
top of this actually worsens the constraints, increasing 
the 95% upper limit on ^m^, from 0.32eV to 0.34eV. 
This is due to the important role the measurement of the 
Hubble parameter plays in the analysis for this model, 
and the tension between these two datasets. The higher 
Hq measurement from Riess et al. [361 gives a smaller 
best fit matter density, and such cosmologies generate 
less structure. As such the damping effect of massive 
neutrinos has to be smaller, and so the neutrino mass 
constraint becomes tighter. With the BAO constraint, 
the matter density is larger, and so more neutrino damp- 
ing allowed, leading to a looser constraint. Combining 
the two datasets leads to a slightly worse fit than might 
be expected because of this tension, and consequently 
does not improve the mass constraint. 

There are two possible effects of neutrinos we could 
consider. Here we have addressed the mass of the neu- 
trino, but recently there has also been a slight preference 
for more than three neutrinos species [iStiSU] . We con- 
sider this subject in an upcoming paper. 



6. Running of the power spectrum 

The primordial power spectrum of density perturba- 
tions may not be a simple power law, where the spectral 
index, n^, is independent of wavenumber, k. Instead the 
power spectrum may have some running [511, such that 
it is now parameterised as 



P{k) = A 



(24) 



where riiun = din rig /din A:, and k^, is some pivot scale. 
The presence of a large running of the spectral index, 
'^run > {iT-s ~ 1)^, is a departure from the predictions 
of ordinary single-field slow-roll inflation, and requires a 
more complex inflationary history (e.g., multiple scalar 
fields 1521 or a non-miminal coupling between the inflaton 
and Ricci tensor [53.). 

We consider a model where the running is allowed to 
vary, rather than being held fixed at Tirun = 0. The 



marginalised constraints are given in table VIII 



The two-dimensional parameter constraint contours 
for {rig, rirun} using CMB data alone, and combining 
CMB with the WiggleZ and other LSS data are shown 
in Fig. [TSj The WiggleZ data drives the running to 
slightly more negative values, though CMB -I- WiggleZ 
is consistent with zero running within two-sigma (nj-un — 
— 0.04 ± 0.022). This result is driven by the large scale 
structure data, as the addition of this extra parameter 
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FIG. 18: (color online). The two-dimensional constraints in 
the {hs, Krun} planc. Adding the WiggleZ data reduces the 
size of the error bars with a best fit value for the running, 
which is more negative than for WMAP alone (rirun = —0.033 
for CMB only, whereas it is n^un = -0.040 for CMB + Wig- 
gleZ). 



allows the combination of CMB and WiggleZ to be con- 
sistent with a slightly larger value of the matter density, 
= 0.332 ± 0.039. This is similar to, but less extreme 
than an effect we have seen before, in section [V D 3[ when 
considering curvature. Combining with other datasets re- 
duces the significance of this result, as the matter density 
is pushed back to lower values. Similarly, measurements 
of the matter power spectrum turnover would also al- 
lowed such large negative running models to be falsified. 



7. Tensors 

We also consider models where the primordial power 
spectrum is composed of both scalar and tensor pertur- 
bations. Tensor perturbations are generated by primor- 
dial gravitational waves from cosmological inflation (for 
reviews see [54H58]). We parameterise the tensor con- 
tribution through the tensor to scalar amplitude ratio, 
defined as 



At 



(25) 



where As is the amplitude of scalar perturbations at some 
pivot scale fc* , and At is the amplitude of tensor pertur- 
bations at that same scale. If infiation is the only mech- 
anism for generating both the scalar and tensor pertur- 
bations (which is the case in single-field, slow-roll infla- 
tion), there is a consistency relation between the tensor 
to scalar ratio and the spectral index of the tensor power 
spectrum, n^, given by 



-SriT ■ 



(26) 
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FIG. 19: (color online). The two-dimensional constraints in 
the {hs, r} plane, for a number of data combinations (see 
legend). We see that that the CMB -I- WiggleZ improves on 
the upper limit from CMB alone, and gives about the same 
limit as CMB -I- LSS. Combining aU three, CMB -I- WiggleZ 
-I- LSS gives the best constraint. 



In this analysis we have assumed the inflationary con- 
sistency relation, and only allowed the tensor to scalar 
ratio parameter to vary, as well as the other cosmological 
parameters in the ACDM model. 



The marginalised constraints are given in table VIII 
The two-dimensional parameter constraint contours for 
{us, r} are shown in Fig. 19 We do not detect the ten- 



sor contribution, but do improve the 95% upper limit on 
it from r < 0.36 (CMB alone) to r < 0.21 (CMB + Wig- 
gleZ). This is very similar to constraints that come from 
CMB combined with other LSS data, and the combined 
CMB -I- WiggleZ -|- LSS data give the best constraints, 
of r < 0.18. 



E. Comparison of all models 



In table |IX[ we give the best fit log- likelihood (£) val- 
ues for the different models and dataset combinations 
considered in this paper. The increase in — log(£) when 
the extra datasets are added to the CMB are consistent 
with what is expected for Gaussian data-points with that 
many degrees of freedom (0.5 for Hq, 1 for BAO, 236 for 
SN-Ia, 406 for WiggleZ). 

While we do not use the log-likelihood as a model selec- 
tion statistic, we notice that the WiggleZ data (in combi- 
nation with the other data we have used) is equally likely 
when considered against all models. The only exception 
is the wGDM model (where the Universe is allowed to 
be nonflat), which provides a much better fit than any 
of the other models when compared against the CMB -f 
WiggleZ data only. However, this improvement in fit is 
removed when another distance data set is added, and so 
is likely due to some undiagnosed systematic or simply 



some random coincidence. 

Since the best fit log-likelihoods of all the models are 
approximately equal for a given data combination, any 
model selection power will come from the prior prob- 
ability. A simple model selection criteria would be to 
approximate the prior by simply penalising each model 
best fit log-likelihood by the number of free parameters 
in the model (e.g., the Akaike Information Criteria, or 
AIC [59]). In this case, the model with the least free pa- 
rameters could be considered the best, which here would 
be flat ACDM. This conclusion is consistent with we see 
from the parameter constraints, where the limiting case 
of flat ACDM is always a good fit. 



VI. CONCLUSIONS 

We have compared the measurement of the galaxy 
power spectrum from the WiggleZ survey to predictions 
for the underlying matter power spectrum to conduct a 
likelihood analysis of a number of cosmological models. 
We have modelled the non-linear structure formation and 
observational effects that go into making these predic- 
tions, choosing between a number of different alternative 
models from the literature and testing against numeri- 
cal simulations specifically created to match the galaxies 
measured by the WiggleZ survey. We found that the best 
performing model was one similar to that used for analy- 
sis of the SDSS LRG sample [55], but where the damping 
of the BAO peak is removed. 

We applied this non-linear approach to measurements 
of the WiggleZ power spectrum, and used data from the 
four redshift bins and seven regions to constrain the cos- 
mological parameters. We performed the analysis using 
the CosmoMC analysis package, with an additional mod- 
ule to compute the likelihood of the WiggleZ P{k). We 
make the module and necessary data products available 
with this paper. 

The results from a number of different possible cosmo- 
logical models are as follows: 

ACDM model: In combination with the CMB, the 
WiggleZ data improves constraints on the matter den- 
sity (ri„i = 0.290 ± 0.016) and amplitude of fluctuations 
(cTg = 0.825 ±0.017). 

Dark energy and Curvature: The CMB -|- WiggleZ 
data combination alone is not powerful at constraining 
the dark energy equation of state {w = —0.53 ± 0.29) or 
the curvature parameter {flk = —0.051 ± 0.028. This is 
because introducing an extra dynamical degree of free- 
dom creates a degeneracy with the matter density, fJ^, 
reducing the constraint on this extra parameter. Dis- 
tance measurements from BAO and SN-Ia provide an 
independent and complimentary measurement that re- 
moves this degeneracy, reducing the errors to flk = 
-0.001 ±0.005 on the curvature and w = -1.008 ±0.085 
on the equation of state (using CMB + WiggleZ + BAO 
+ Hq data). 

Massive neutrinos: The CMB + WiggleZ data pro- 
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Model Number of CMB only CMB + WiggleZ + Ho + BAO + Ho+ BAO + SN-Ia + LSS 

Parameters 



Best fit — log(£) values 



Flat ACDM 


7 


3734.1 


4140.9 


4141.8 


4141.6 


4142.6 


4353.2 


4164.4 


Flat wCDM 


8 


3734.0 


4140.3 


4141.6 


4141.6 


4141.9 


4352.8 




ACDM 


8 


3733.9 


4140.6 


4141.5 


4141.5 


4142.6 


4353.2 




wCDM 


9 


3734.2 


4136.5 


4141.4 


4141.4 


4142.5 


4351.4 




Flat ACDM + 


8 


3734.0 


4140.8 


4141.9 


4141.5 


4142.6 




4164.0 


Flat ACDM + n^un 


8 


3733.6 


4139.6 


4141.1 


4140.5 


4141.8 




4162.9 


Flat ACDM + r 


8 


3734.0 


4140.9 


4141.9 


4141.7 


4142.7 




4164.4 


log(£) differences (to Flat 


ACDM) 
















Flat ACDM 


7 


0.00 


0.00 


0.00 


0.00 


0.00 


0.0 


0.0 


Flat wCDM 


8 


-0.10 


-0.62 


0.01 


-0.04 


-0.67 


-0.48 




ACDM 


8 


-0.16 


-0.33 


0.11 


-0.07 


0.05 


-0.04 




wCDM 


9 


0.15 


-4.41 


0.14 


-0.20 


-0.03 


-1.86 




Flat ACDM + E 


8 


-0.06 


-0.11 


0.12 


-0.16 


0.06 




-0.34 


Flat ACDM + n,un 


8 


-0.48 


-1.30 


-0.75 


-1.06 


-0.74 




-1.43 


Flat ACDM + r 


8 


-0.05 


0.04 


0.07 


0.05 


0.10 




0.05 



TABLE IX: Best fit — log{£) values for the different cosmological models and datasct combinations, and the differences between 
those best fit — \og(£) values and the best fit Flat ACDM — log{£). For the second half of the table, a negative value represents 
a better fit of the data to the model, and a positive value a worse fit. 



vides a powerful constraint on the neutrino mass, giving 
an upper limit on the neutrino mass (summed over the 
three species) of 0.58eV (95% confidence). However this 
result is still somewhat degenerate with the matter den- 
sity. Other distance measurements reduce this degen- 
eracy, giving our best upper limit of 0.32 eV (CMB -|- 
WiggleZ + Ho). 

Primordial Power spectrum and gravitational waves: 
The CMB -|- WiggleZ data does not provide any bet- 
ter constraint on the spectral index than is provided by 
the WMAP data by itself. However, by constraining the 
matter density, we can get better constraints on param- 
eters such as the amplitude of fluctuations, erg. By com- 
bining with other large-scale structure data, we can also 
improve measurements on the running of the spectral in- 
dex nrun = -0.033 ± 0.018 (CMB -|- WiggleZ + LSS) or 
the tensor to scalar amplitude ratio r < 0.18 (CMB -|- 
WiggleZ + LSS). 

To summarise, all our results were consistent with the 
standard flat ACDM cosmological model, though we have 
significantly improved constraints on a number of param- 
eters. 

We see that many of the models created degeneracies 
with rim that the CMB -I- WiggleZ data combination can- 
not break, giving multiple or extended likelihood peaks, 
and best fit parameter values with high matter densities 
and low values of Hq. This is because the power spec- 
trum cannot distinguish models with an extra degree of 
freedom if f2m takes a larger value (for example, curva- 
ture, massive neutrinos and running of spectral index). 
These high matter density peaks can be removed if ex- 
tra data is added, but could also be falsified if the power 



spectrum were able to accurate measure the turnover at 
large scales. The WiggleZ data has been used to deter- 
mine the turnover position and so measure the matter 
density f2„i in a separate publication [44 . 

The analysis presented in this paper has demonstrated 
the importance of modelling non-linear structure forma- 
tion, bias, and redshift space distortions. The importance 
of this modelling will only increase as future galaxy sur- 
veys probe the matter power spectrum at higher precision 
and down to smaller scales. 
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Appendix A: Fitting fomulae for GiggleZ Power 
Spectum 

We fit the simulated power spectra from GiggleZ using 
a fifth-order polynomial in k (in units of /i~^Mpc), for 
the logarithm of the power spectrum, P(fc) (in units of 
h'~'^ Mpc'^), in each of the four redshift bins. In order 
to generate mock WiggleZ catalogues from the GiggleZ 
simulations at a particular redshift, we selected the range 
of halo groupings with large-scale clustering bias closest 
to the WiggleZ sample under analysis, and applied the 
WiggleZ selection function to the mock dataset. 



2 


^0 


Al A2 


A3 


A4 


^5 


0.22 4.619d0 


-13.7787 58.94 


-175.24 


284.321 


-187.284 


0.41 


4.63079 


-12.6293 42.9265 


-91.8068 


97.808 


-37.633 


0.60 


4.69659 


-12.7287 42.5681 


-89.5578 


96.664 


-41.2564 


0.78 


4.6849 


-13.4747 53.7172 


-145.832 


216.638 


-132.782 



TABLE X: Coefficients for the polynomial fitting formula for 
the GiggleZ power spectrum. 



The resulting fitting formula is: 

logio(P(A:)) = Ao + Aik + A2k^ (Al) 
+A3k^ + Aik"^ + A^k^ 

with different coefficients for each of the four redshift 
bins. The coefficients are given in table |Xj and built-in 
in the CosmoMC module. 



